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Abstract. In the field of radial basis functions mathematicians have been 
endeavouring to find infinitely differentiable and compactly supported radial 
functions. This kind of functions is extremely important. One of the reasons 
is that its error bound will converge very fast. However there is hitherto no 
such function which can be expressed in a simple form. This is a famous 
question. The purpose of this paper is to answer this question. 
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1 Introduction 

In this paper we raise a new function 

:= { o if in > i W 

where a > 0, x G R n and ||a;|| denotes the Euclidean norm of x. The infinite 
differentiability of $ can be shown by routine investigation, and we omit it. 
What's difficult is to show its positive definiteness. In order to make it useful 
in scattered data interpolation, this problem has to be overcome. This is the 
main theme of the next section. 
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2 Positive Definiteness 



It's well known that positive definiteness is just conditional positive definite- 
ness of order zero. Here we adopt the definition of pQ for conditional positive 
definiteness. Before showing that the map defined in (1) is positive definite, 
let's define 

In other words, 0(||x||) = for x G R n , where a > is a constant. In 

order to make this paper easier to understand, let's analyze the derivatives 
of 0(V?) of order up to four. 

It's easily seen that with u = tan^t, 

. = e -«(i+«) 2 (- a )|(l +M 2)[ 2n + 2] 

• J^(Vt) =e- Q ( 1+ ") 2 «(f) 2 (l + M 2 )[4a(l + u) 2 (l + M 2 )-2(l + n 2 )-4(l + 

• = e- Q ( 1+ ") 2 (-«)(f ) 3 (1 + w 2 )[8a 2 (l + w) 3 (l + w 2 ) 2 - 12a(l + 
w)(l + it 2 ) 2 - 24cra(l + w) 2 (l + w 2 ) + 8(1 + u 2 )u + 4(1 + 2w)(l + u 2 ) + 
8(l + w)w 2 ] 

• ^(p(Vt) = e- a ( 1+ ") 2 «(f) 4 (l + M 2 )[16a 3 (l + u)\l + m 2 ) 3 - 48a 2 (l + 
w) 2 (l + it 2 ) 3 - 96a 2 (l + m) 3 m(1 + w 2 ) 2 + 136a(l + u)u{l + u 2 ) 2 + 8a(l + 
u)(l+2M)(l+M 2 ) 2 +64a(l+n)V(l+u 2 ) + 12a(l+M 2 ) 3 +24a(l+u) 2 (l + 
3u 2 )(l + u 2 ) - 104u 2 (l + u 2 ) - 32m(1 + u 2 ) - 16(1 + u 2 ) - 16u 3 (l + u)]. 

By induction, we find 

^(Vl) = e-^ +u y\~iya^y(l + u 2 )F J (a,u). (3) 
where Fj(a,u) is a polynomial of a, u with integer coefficients. Note that 

j-2 3j-3 

Fj-(a,u) = o j a J '- 1 u 3 t J "- 1 5 +1 + £ £ a w aV 

fc=0 «=0 

where a., G Z + and a^i G Z for all k,l. 

Lemma 2.1 Suppose Fj is defined as in (3) forj > 1. There exists a number 
atj > such that Fj(a, u) > for all a > ctj and u > 0. 
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Proof. Obviously Um 



a 



oo 



Fj(a, u) 



oo since aj > 0. For any M > 0, 



1 1 



oo 



there exist um > and olm > such that Fj(a,u) > M whenever u > um 
and a > %. Let M = 0. Then Fj(a,u) > for all u > uq and a > a7o. 
For each w G [0,wo], Fj(a,u) is a polynomial of a. Let it be p u (oc). Since 
lirria^ooPu^a) = oo, there exists a u > such that > whenever 

a > a M . Let g(u) be the infimum of such ot u . Then g(u) is obviously a 
continuous function on [0,Mq]. Let (3j = sup u£ [ ^ uo ]g(u). Then, for all a > 
(3j, Fj(a,u) > if u G [0, it ]. Now, let acj = max{a , {3j}. Then Fj(a,u) > 
for all a > otj and u > 0. (j 

Theorem 2.2 For any positive integer n, there is a positive number M n such 
that for all a > M n , the function $ defined in (1) is strictly positive definite 
on R d for 1 < d < n. 

Proof. By Theoreml of pQ and Theorem2.2 of |2], for each 1 < d < n, if 
(-l) J ^0(v / t) > Ofor j = 0,1,. . .,/on (0,oo), where / = [d/2] + 2, then 
is strictly positive definite on R d . Let nid = max{a\, . . . , a{\ where oti, . . . , ai 
are defined in Lemma2.1. We have (— l)- 7 J^-0(-y/t) > for j = 0, 1, . . . , / on 
(0, oo) if a > nid. Now, let M n = max{mi, . . . , m n }. The theorem then fol- 
lows immediately. (j 

Example. The function defined in (1) is strictly positive definite on R d 
for 1 < d < 4 whenever a > 2. 



Proof. By Theoreml of pQ and Theorem2.2 of [2], it suffices to show that 
§?4>(Vi) > for j = 0, 1, . . . , 4 on (0, oo), where is as in (2). For 
j = 0, 1, 2, it's trivial. We begin with j = 3. 
For j = 3, 



F 3 (a,u) = 8a 2 (l + w) 3 (l + w 2 ) 2 -12a(l + w)(l + w 2 ) 2 

-24cm(l + u) 2 (l + u 2 ) + 8(1 + u 2 )u + 4(1 + 2u)(l + m 2 ) 
+8(1 + u)u 2 




a(-f(l + u 2 )F 3 (a,u) 



where u = tan^t and 
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= 4a(l + u)(l + M 2 )[2a(l + m) 2 (1 + m 2 ) - 3(1 + m 2 ) - 6«(1 + u)] 
+8(1 + m 2 )m + 4(1 + 2u)(l + m 2 ) + 8(1 + u)u 2 . 

The terms in the brackets equal 

2a(l + 2m + m 2 )(1 + m 2 ) - 3(1 + m 2 ) - 6u(l + «) 
= 2a(l + m 2 ) + 2aM 2 (l + m 2 ) - 3(1 + m 2 ) + 4au(l + u 2 ) - 6u(l + «) 
> 4cm(l + m 2 ) - 6u(l +«) + 2aM 2 (l + m 2 ) 
= 4au + 4au s — 6m — 6u 2 + 2au 2 + 2au 4 
= 2u[2au 2 + (a - 3)m + 2a - 3] + 2aM 4 

For 3 < a, (a — 3)m + 2a — 3 > 0. For 2 < a < 3, and m > 1, it's obvious 
that 2cm 2 + (a — 3)m + 2a — 3 > 0. For 2 < a < 3, and < u < 1, we have 
| (a — 3)m| < 1 and 2a — 3 > 1. This leads to 2cm 2 + (a — 3)m + 2a — 3 > 0. 
In any case, F 3 (a, u) > for < u < oo and a > 2. 
We now investigate the case j = 4. For j = 4, 

(-l) 4 ^0(v^) = e-^) 2 a(|) 4 (l + M 2 )F 4 (a, M ) 
= e -°a+«) 2 a (J) 4 (l + M 2 ) 2 

where 

= 16a 3 (l+M) 4 (l+M 2 ) 2 -48a 2 (l+M) 2 (l+M 2 ) 2 

1 + u 2 

-96a 2 (l + ufu(l + u 2 ) + 136a(l + u)u{\ + m 2 ) 

+8a(l + u)(l + 2u)(l + m 2 ) + 64a(l + m) 2 m 2 

+ 12a(l + m 2 ) 2 + 24a(l + m) 2 (1 + 3m 2 ) - 104m 2 - 32m - 16 

-i 6 ." 3(1 + "l 

1 +M 2 

Let G(a,u) := F *+^ ■ The crucial part of G(a,u) is the first three terms. 
The three terms equal 

8a 2 (l + m) 2 (1 + m 2 ){(1 + M 2 )[2a(l + m) 2 - 6] - 12(1 + m)m}. 

Let g a (u) = (l+M 2 )[2a(l+M) 2 -6]-12(l+M)M. It's easily seen that g a {2) > 
for a > 2. Now, g' a (u) = 2u[2a(l +m) 2 - 6] + (1 + m 2 ) [4a (1 +u)] - 12(2m+ 1). 



F 4 (a, u) 

1 + M 2 
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For a = 2, g' 2 (u) = 8(1 + u)(l + m 2 ) + 8(1 + m) 2 m - 36m - 12 > if u > 2. 
Obviously g' a (u) > for a > 2 and u >2. These lead to <? q (m) > for a > 2 
and u > 2. It's then easy to see that G(a, m) > for a > 2 and u > 2. 

Before analyzing G(a,u) for a > 2 and < u < 2, let's find J^G(o:,m) 
and -^G(a,u) first. 

-j-G(cx, u) = 16a 3 4(l+M) 3 (l + M 2 ) 2 + 16a 3 (l + M) 4 2(l + M 2 )2M 



du 



-48« 2 2(1 + u)(l + m 2 ) 2 - 48a 2 (l + m) 2 2(1 + u 2 )2u 

-96a 2 3(l + m) 2 (m + m 3 ) - 96a 2 (l + m) 3 (1 + 3m 2 ) 

+ 136cm(l + u 2 ) + 136a(l + u)(l + 3m 2 ) 

+8«[8m 3 + 9m 2 + 6m + 3] + 128a(l + m)m 2 + 64a(l + m) 2 2m 

+24a(l + m 2 )2m + 24a2(l + m)(1 + 3m 2 ) + 24a(l + m) 2 6m 

-208, - 32 - 16 ■ (1 + »')(3 M ' + 4 U ') -~'(1 + M )2u 

(1 + M ) 

64a 3 (l + m) 3 (1 + m 2 ) 2 + 64a 3 (l + m) 4 (1 + m 2 )m 
-96a 2 (l + m)(1 + m 2 ) 2 - 480a 2 (l + m) 2 (m + m 3 ) 
-96a 2 (l + m) 3 (1 + 3m 2 ) + 184au(l + m 2 ) 
+184a(l + m)(1 + 3m 2 ) + 8«[8m 3 + 9m 2 + 6m + 3] 
+128a(l + m)m 2 + 272a(l + m) 2 m - 208m 
_o 9 _ 1fi (1+m 2 )(3m 2 + 4m 3 )-m 4 2(1 + m) 

62 (1 + M 2 ) 2 W 



—G(a,u) = 48« 2 (1 + m) 4 (1 + m 2 ) 2 - 96a(l + m) 2 (1 + m 2 ) 2 
da 

-192a(l + m) 3 m(1 + m 2 ) + 136(1 + u)u(l + u 2 ) 

+8(1 + m)(1 + 2m)(1 + m 2 ) + 64(1 + m) 2 m 2 

+ 12(1 + m 2 ) 2 + 24(1 + m) 2 (1 + 3m 2 ) (5) 

The first two terms of -£^G(a,u) equal 

48a 2 (l + 4m + 6m 2 + 4m 3 + m 4 )(1 + m 2 ) 2 - 96a(l + m) 2 (1 + m 2 ) 2 
= 48a 2 (l + 2m + m 2 )(1 + m 2 ) 2 + 48« 2 (2m + 5m 2 + 4m 3 + m 4 )(1 + m 2 ) 2 

-96«(1 + m) 2 (1 +m 2 ) 2 
> 48a 2 (2m + 5m 2 + 4m 3 + m 4 )(1 + m 2 ) 2 if a > 2 and u G (0, oo). 
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Hence 

-j-G(a, u) > 48a 2 (2w + 5w 2 + 4w 3 + w 4 )(l + w 2 ) 2 - 192a(l + w) 3 w(l + w 2 ) 
da 

+136(1 + u)u(l + u 2 ) + 8(1 + u)(l + 2u)(l + m 2 ) 

+64(1 + u) V + 12(1 + m 2 ) 2 + 24(1 + w) 2 (l + 3m 2 ) (6) 

The first four terms of (6) equal 

(1 + w 2 )8[6a 2 (2w + 5m 2 + 4m 3 + w 4 )(l + u 2 ) - 24a(l + ufu 
+17(l + u)u+(l + u)(l + 2u)] 
= (1 + w 2 )8[6aV + 24a V + (36a 2 - 24a)w 4 + (36a 2 - 72a)w 3 
+ (30a 2 - 72a + 19)w 2 + (12a 2 - 24a + 20)w + 1](> if u > 1) 

> (l+w 2 )8[(30a 2 -72a + 19)w 2 + 20w + l] 

> (1 +u 2 )8[-5u 2 + 20u + 1] 

> (1 +« 2 )8[-5w + 20w + 1] i/ < u < 1 

> (1 +m 2 )8[15m + 1] 

> i/ < u < 1 (7) 

Thus 

— G(a, u) > if a > 2 and w G (0, oo). (8) 
da 

Now, let a = 2. Then 

^G(2,u) = 2\l + uf(l + u 2 ) 2 + 2\l + u) A (l + u 2 )u 

-3 • 2 7 (1 + u)(l + m 2 ) 2 - 480 • 4(1 + m)V 1 + u 2 ) 
-(2 8 + 2 7 )(1 + uf(l + 3m 2 ) + u(l + u 2 ) [2 7 + 3 • 2 6 + 2 5 + 2 4 ] 
+ (1 + u)(l + 3« 2 )[2 7 + 3 • 2 6 + 2 5 + 2 4 ] 
+ 16(3 + 4«)(1 + u 2 ) + 32(1 + 3u + 2w 2 )w + (1 + w)w 2 • 2 8 
+ (1 + w) 2 w(2 8 + 3 • 2 6 + 3 • 2 5 ) - 208m - 32 
(1 + u 2 )(3u 2 + 4m 3 ) - u 3 (l + u)2u 
{l + u 2 y 



~ 16 7V— ( 9 ) 



We are going to analyze the value of G(2,u) for 1 < u < 2. Note that for 
1 < u < 2, (1 +u) A > 4(1 +u) 2 . By this the first five terms of (9) can easily 
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be analyzed as follows. 

2 9 (1 + m) 4 (1 + u 2 )u - 480 • 4 • (1 + m) 2 m(1 + m 2 ) 
> 4 • 2 9 (1 + m) 2 (1 + u 2 )u - 1920(1 + u) 2 u(l + u 2 ) 
= 128(1 + m) 2 (1 + u 2 )u (u G [1,2]) (10) 

Then 

128(1 + m) 2 (1 + u 2 )u - 3 • 2 7 (1 + u)(l + m 2 ) 2 
= 128[(1 + m 2 ) 2 m + 2(1 + u 2 )u 2 } - 3 • 2 7 (1 + u)(l + u 2 ) 2 
= -256(1 + m 2 ) 2 m- 384(1 + u 2 ) 2 + 256(1 + m 2 )m 2 . (11) 

Now, 

2 9 (1 + m) 3 (1 + u 2 ) 2 - (2 8 + 2 7 )(1 + m) 3 (1 + 3m 2 ) 
= 2 9 (1 + uf{\ + m 2 ) 2 - (2 8 + 2 8 )(1 + u)\l + 3m 2 ) + 2 7 (1 + uf{\ + 3m 2 ) 
= 2 9 (1 + m) 3 [1 + 2m 2 + u A - 1 - 3m 2 ] + 2 7 (1 + m) 3 (1 + 3m 2 ) 
= 2 9 (1 + m) 3 [m 4 - m 2 ] + 2 7 (1 + m) 3 (1 + 3m 2 ) 

> 2 7 (1 + m) 3 (1 + 3m 2 ). (m G [1,2]) (12) 

The sum of (11) and (12) is 128m 5 + 1024m 4 + 256m 2 + 128m - 256 which is 
larger than if u G [1,2]. Thus, the sum of the first five terms of (9) is 
positive. It's trivial to see that the sum of the other terms is also positive. 
Hence 

|G(2,»)>0 (13) 

for m G [1,2]. Since G(2, 1) > 0, by (8) and (13), we get G(a,u) > for 
a > 2 and u G [1,2]. 

The next task is to analyze G(a, u) for a > 2 and u G (0, 1). We first let 
a = 2 and u G (0, 1). In this case, 

G(2,u) = 16-2 3 (1+m) 4 (1+m 2 ) 2 - 48-2 2 (1 +m) 2 (1+m 2 ) 2 

-96 • 4(1 + m) 3 (m + m 3 ) + 136 • 2(m + m 2 )(1 + m 2 ) 

+ 16(1 + m)(1 + 2m)(1 + m 2 ) + 128(1 + m) 2 m 2 

+24(1 + m 2 ) 2 + 48(1 + m) 2 (1 + 3m 2 ) - 104m 2 

m 3 (1 +u) 
-32m - 16 - 16 • V V 

1 + M 
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= 2 V + 2V + (5 • 2 6 + 2 7 )m 6 + 280m 4 + 96m 3 - 120m 2 




> 8m + 8 - 16(m 3 + m 4 ) + 96m 3 + 280m 4 (since u E [0, 1]) 
= 8m + 8 + 80m 3 + 264m 4 



> 8. 



This result together with (8) shows that G(a, u) > for a > 2 and u £ (0, 1). 



Remark. For dimensions higher than four, the value of a can be analyzed 
analogously. It's a bit more tedious, but not difficult. For each increase of 
the dimension by two, the order of the derivative of 4>(Vi) has to be increased 
by one. 
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Our conclusion thus follows. 
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